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DETERMINATION OF LATERAL-STABILITY DERIVATIVES AND TRANSFER-FUNCTION 
COEFFICIENTS FROM FREQUENCY-RESPONSE DATA FOR LATERAL MOTIONS ^ 

By Jambs J. Doneoan, Samuel W. Robinson, Jb., and Ordwat B. Gates, Jk. 


SUMMARY 

A method is 'presented for determining the lateral-stahility 
derivatives, transfer-function coefficients, and the modes for 
lateral motion from frequency-response data for a rigid aircraft. 
The method is based on the application of the vector technique 
to the equations of lateral motion, so that the three equations of 
lateral motion can he separated into six equations. The method 
of least squares is then applied to the data for each of these 
equations to yield the coefficients of the equations of lateral 
motion fi'om 'which the lateral-stahility derivatives and lateral- 
motion transfer-function coefficients are computed. Two n'u- 
merical examples are given to demonstrate the use of the method. 

INTRODUCTION 

In the reduction and generalization of flight-test data, 
whether for loads, stability, or control purposes, the airplane 
stability derivatives and the coefficients of the transfer 
functions are often required'. A great deal of emphasis, 
therefore, has been placed on the development of analytical 
methods for reducing flight data to obtain these basic deriva- 
tives and coefficients. 

A number of recent methods, for example references 1 to 
4, are now available for analyzing longitudinal maneuvers 
and determining the longitudinal-stability derivatives a,nd 
transfer-function coefficients from flight data. References 
1 and 2 present methods of determining the longitudinal- 
stability derivatives and transfer functions dkectly from 
transient data. Reference 3 reduces data for longitudinal 
motion deteimined from the forced-oscUlation technique by 
means of circle diagrams to longitudinal-stabUity derivatives 
and frequency response. Mueller, in reference 5, was one 
of the earliest to use vector representation in the equations 
of longitudinal motion to represent the derivatives and in- 
tegrals of the variables. Schumacher, reference 4, repre- 
sented the frequency responses to longitudinal motion as 
vectors and substituted them into the equations of longi- 
tudinal motion and the transfer functions. He then applied 
the method of least squares to these vector equations, a 
method which he found very effective in determini n g certain 
of the longitudinal-stability derivatives and transfer-function 
coefficients. 

The problem of analyzing lateral motions, however, has 
not received the same amount of attention as that for longi- 
tudinal motion, perhaps because it is more complicated. 
Several analytical investigations have been undertaken and a 


few methods have been proposed such as the circle-diagram 
method (ref. 6), the step-function-response method (ref. 6) 
in which the response of an airplane to a step deflection of 
the rudder or aileron is analyzed, and the free-oscfllation 
method (ref. 7) in which the period and damping of the 
free vibrations of the aircraft due to a pulse-type input are 
analyzed. Since the usefulness of these methods is limited 
by the number of derivatives which can be extracted, there 
is still a need for a more general method of analysis that will 
extract all the significant lateral-stability derivatives from 
flight data. 

It is the purpose of this report to present a method for 
determining the lateral-stability derivatives of a rigid air- 
plane and to illustrate its use by applying it to two examples. 
The method is based on the vector representation of the 
frequency responses to lateral motions. This vector ap- 
proach permits separation of each of the equations of lateral 
motion into a real and imaginary equation. A least-squares 
method is then applied to the data in each of these equations 
or combinations thereof to yield the coefficients of the equa- 
tions of latei’al motion. The lateral-stability derivatives 
and transfer-function coefficients are then determined from 
these coefficients and the known aircraft mass parameters. 

The method is applied to two speciflc examples, one in 
w'-hich the frequency, responses to a rudder input are known 
and one in which the transient responses to aileron deflection 
are known. In the latter case, the frequency responses were 
obtained from the transient motions by two methods and 
the stability derivatives computed.. 

An attempt has been made to schedule the procedure so as 
to reduce the dependence of the results obtained from this 
method on the derivatives that can be least accurately 
obtained from the particular data being analyzed; however, 
fm’ther improvements may be made as further experience is 
gained in the application of the method. 

SYMBOLS 

ay lateral acceleration, ft/sec^ 

Aff, A^,, A,j,, parameters defined by equations (8) 

B^,, Ba^ parameters defined by equations (9) 

h wing span, ft 

(7i lift coefficient, L/qS 



C„ yawing-moment coefficient, N/qSh 


‘ Supersedes NACA TN 3083, 1954. 
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Ct 

Or 

Cl, C-2, C3, 


D 


Fi, F„ F3 

Ci, 63 

s 

Ix 

Iz 

Ixz 


i=V~ 1 

Ki, K3, K3, 


L 

L' 

M 

m 

N 

p or ^ 


r or ^ 


rolling-moment coefficient, L'lqSb 
lateral-force coefficient, Lateral iovc&jqS 
coefficient of transfer functions (defined 
in table II) 

df 

differential operator, ^ 

natural logarithmic base 
forcing-function coefficients representing 
rudder effectiveness (defined in table I) 
forcing-function coefficients representing 
aileron effectiveness (defined in table I) 
acceleration due to gravity, ft/sec* 
moment of inertia about stability A^-axis, 
mkx^, slug-ft^ 

moment of inertia about stability Z-axis, 
mkz^, slug-ft* 

product of inertia refeired to stability 
axes (negative when the positive direc- 
tion of the X principal axis is inchned 
above the flight path, i. e., when ij is 
positive) 

stability-derivative coefficients of the 
equations of lateral motion (defined in 
table I) 

radius of gyration about stability A'’-axis, 
ft 

radius of gyration about stability Z-axis, 
ft 

nondimensionai radius of gyration about 
longitudinal stability axis, 



nondimensionai radius of gyration about 
vertical stability axis. 



nondimensionai product-of-inertia param- 
eter, — 

radius of gyration about principal longi- 
tudinal axis, ft 

radius of gyration about principal vertical 
axis, ft 
lift, lb 

rolling moment, ft-lb 
Mach number 
aii'craft mass, W(g, slugs 
yawing moment, ft-lb 
incremental rolling angular velocity about 
A'-axis, radians/sec 

dynamic pressure, i pV^, Ib/sq ft 

incremental yawing angular velocity 
about Z-axis, radians/sec 
amplitude ratio 
wing area, sq ft 
Laplace transform variable 
time, sec 

incremental component of velocity 1^ 
along the F-axis, ft/sec 


V 

W 

X, Y, Z 


Bai- notation: 


true airspeed, ft/sec 
aircraft weight, lb 
airplane stability axes (see fig. 1) 
angle of attack, radians 
angle of sidelsip, v/V, radians 
aileron control deflection, radians 
rudder control deflection, j’adians 
inclination of principal longitudinal axis 
of inertia with respect to flight path 
(positive when the positive direction 
of the X principal axis is inclined above 
the flight path) 

nondimensionai mass parameter used foi* 
lateral equations, m/pSb 
air density, slugs/cu ft 
time parameter, mIpSV, sec 
phase angle, radians 
angle of roll, radians 
angle of yaw, radians 
angular frequency, radians/scc 

bar over letter represents maximum 
value 

bars on sides of symbol represent ab- 
solute value 


Matrbc notation: 

[ ] square matrix 

{ } column matrix 

The lateral-stability derivatives are expressed by sub- 
script notation as, for example: 


b/S ’ pb’ ^ 1-6 ■ 

Phase angles are also indicated by subscript notation as 
is the phase angle between the input rudder deflec- 
tion and the output sideslip angle. 

DEVELOPMENT OF METHOD 
EQUATIONS OF LATERAL MOTION 

The equations of lateral motion based on the stability axes 
(as shown in fig. 1) and on the usual assumptions of linearity, 
small angles, and maneuvers which start from a level-flight 
condition are usually wiitten in the following form: 

(a) For the sideslipping motion, 

qD+Ki)^-K3<!>+D4>=^FMt) ( 1 ) 

or, since a„=y(D/3-h£)^— isTai^)), 


■yaf-^Ki^=Fi5t{f) (la) 

(b) For the rolling motion, 

^ 3 / 3 + {KJD^^K3m=FMt) + 6-Mt) ( 2 ) 

(c) For the yawung motion, 

-m- {K3D^^KM+ {D^^Kim^F3b,{t) + CM) ( 3 ) 
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Figure 1. — Stability axis system showing positive directions. 


where the K, F, and G coefficients are defined in table I. 
The forcing functions of these equations are written for rud- 
der inputs and aileron inputs combined; for the sake of sim- 
plicity in developing and demonstrating the method, only 
rudder inputs Fidr(t), F^Srit), and FiSrit) are considered. 
The method, however, applies equallj’- well to both aileron 
and rudder inputs as is shown subsequently. 

VECTOR INTERPRETATION OF FREQUENCY RESPONSE 

In manj' current studies of au’plane dynamics, the fre- 
quency response of the aircraft is determined from flight data. 
Having the data available in the frequency plane offers cer- 
tain advantages over having the data in the time plane, since 
use can then be made of the vector interpretation of fre- 
quency 1 ‘esponse. The development which follows will as- 
sume that the data to be analyzed ai’e available in frequency- 
response form. In the event that only transient-response 
data are available, the transformation to the frequency plane 
can be made from a selection of one of the several methods 
compared in reference 8. 

The vector technique is applied to the determination of 
lateral-stability derivatives in the followmg manner. If an 
input 5, to a linear system (a system described by a linear 
differential equation) is considered to have a sinusoidal vari- 


TABLE I.— DEFINITIONS OF THE COEFFICIENTS OP THE 
EQUATIONS OF LATERAL MOTION OF AIRPLANE 


Stability-derivative coefflcients : 

Control-effectiveness coefficients: 

(-s) 

K,=Cl Y^=glV 
Kz-Ci^ [ 2(fci/6)M] 

4r(lbW6)J 

K.^IxzHx 

4r(fcx/6)^ 

TT r 

' 2(fcz/6)M 

Ka= Ixzjlz 

4r(fcz/6)“ 

ET n 

'»r2(fcx/6)M 

n n M6 

2(fcz/&)M 


imaginary oxis 



(a) I 

I 

(a) Vector diagram showing relative magnitudes and positions of input 
Sr and output /3, and their derivatives. w=0.7 radian/sec. 

Figure 2. — Vector representation of lateral variables. 

ation of frequency w, it can be considered to be a vector of 
unit magnitude lying along the real axis of the complex plane 
as shown in figure 2(a). The first derivative of the sinus- 
oidal input DSr is then obtained by multiplying the amplitude 
by w and rotating the resulting vector 90° counterclockwise 
in the complex plane which is equivalent to advancing the 
phase angle by 90°. Each succeeding higher derivative is 
found by multiplying by w and rotating the vector 90° in 
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the complex plane; thus, 


fi,=l 

(4) 

Z?S,=l’w 

(5) 

Zl-5,= — w* 

(6) 


The steadj'-state output or response of the linear s3’-stem 
to a sinusoidal input may be considered to be a vector of mag- 
nitude equal to an amplitude ratio R and having a dhection 
angularty displaced from the real axis by a phase angle 
Thus, in the complex plane, the vector for sideslip may be 
represented as 




1 

Sr 

e = 

Sr 


^cos sin 




(7) 


i 

! Imoginary axis 



By definition. 


Aff= 

1 

COS 

A,j,= 

± 

■* 

COS 

«r 

> 

A^= 

1 

Sr 

cos$,^^ 

III 

at 

.4 

COS 

B,^ 

£ 

Sr 

sin 


4 

sin 







74= 

1 

Sr 



5: 

4 

sin \ 


The derivatives of jS are represented 'as 


Z?j8=w 


e 


sin +i cos 


and 






= W“ 


(-cos sin #,3^ 


= — (o^Ap — iorB^ 


(81 


(9) 


( 10 ) 


( 11 ) 


Similarly, for loll (see fig. 2(b)), the equations corresponding 
to equations (7), (10), and (11) would be 


for yaw, 


<li=Aifi-\-iB^ 

(12) 

D<f>=—o3B^A'i(^-i^4> 

(13) 

D^<l> = — (d^A^ — iu'~B^ 

(14) 

A^-\~iB^ 

(15) 

D\p= — caB^Ai(^A^ 

(16) 

Z?V— — (li^Af — %(AB^ 

(17) 


( 

I 

(b) Vector diagram showing relative magnitudes and positions of inpul 
Sr and response vectors (p, and 4/. 

Figure 2. — Concluded, 
and, for lateral acceleration, 

0.il=Aay-\-iBa^ ( 18 ) 

Substituting equations (4) to (18) into equations (1), (la), 
(2), and (3) and equating the real and imaginary values in 
each equation yields the following eight simplified equations: 
From equation (1): 

ii:i5^-/v2g*=-co(^l3-h.-4) (10) 

and 

( 20 ) 

From equation (la): 

VK:B^=-Ba^ ( 21 ) 

T7vV4=rF,-X (22) 

From equation (2): 

K,B^+KiU)A^+K,u^B^-Iu<»A^=o}-B^ (23) 

and 

KiAff-KiwB^+Ks<^^A^+K,^Bt-F2=u-A^ (24) 

From equation (3) : 

(2.'5) 

and 

— KiAff Agto44 Kg(i)B,p — Aiflw/4 — Fz = uv 4 (20) 

In equations (19) to (26) the and 74,^,# ,o^ terms 

will be available at the particular values of w from the fre- 
quenej'-response curves. It will be assumed that the A'j 
term is determined from the velocitj' measuiement and the 
Ks and Ks terms, which are equal to the ratios Ixzllx and 
Ixztlz, respectively, are known from eitlier weight and 
balance calculations or measurements so that only the values 
of Ai ,3,4,6,7.9.10 and ^1,2.3 are to be determined from the data. 
(It should be noted that even though /ij and As are assumed 
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to be known a need still exists for a simple method of deter- 
mining aircraft moments of inertia and the location of the 
principal longitudinal axis.) 

Equations (19) to (26) are now used to compute the K 
and F coefficients. The advantage of separating each of 
the equations of lateral motion into real and imaginary 
equations now becomes apparent. The two derivatives 
most difficult to determine are Ci^ and contained in the 
Ki and Kg coefficients, respectively. The Kg coefficient may 
be eliminated between equations (23) and (24); similarly 
the Kg coefficient may be eliminated between equations (25) 
and (26). The equations thus obtained, when used to 
compute the unknown coefficients, result in better condi- 
tioned matrices and more accurate results. In the usual 
case, it has been found unnecessary to perform the elimina- 
tion of the Kg coefficient. 

As a first step in solving for the K and F coefficients, 
equations (19) to (26) are fitted to the frequency-response 
data over a range of values of w by the method of least 
squares. The theory of least squares is well known and is 


derived in many textbooks (e. g., ref. 9) and no attempt is 
made here to repeat its derivation. 

The application of the least-squares method to equations 
(19) and (26) converts these equations to their computational 
form. It is recommended that, if the lateral-acceleration 
frequency response is available, equations (21) and (22) be 
used in preference to equations (19) and (20) for the com- 
putation of Ki and Fi. 

Equation (19) may be expressed simply as 

(tin 

K,=^^^ (27) 

S m)if 

and yields Ki. Similarly, equation (20) can be expressed as 

1 Q>n 

Ei=- y) {K^A^-K2A^-o,B^-<cB^)i (28) 

and yields Ej. 


For the same variables K and Fj, the least-squares forms of equations (21) and (22) for use with the lateral-acceleration 
frequency response are: 


■^1— 6., 

S [(yB,)A^ 


(29) 


Equation (25) can be expressed as 


i)[(-5,),]^ 

<Jn * 

;=3<oi 

2 {<^Ai)j{—B^j 


V. 


' f^ibi^B^-Kgbi^BM-B?)! 

}-m 

<j3n 

j—ia\ 

i) [{-b,A,)jf 

Wn 

-< 

Kg 



S {~B^j{uA,i)j 

_ j—<a\ 

<«>i» 

(— wA^)j(wA^)^ 



..Ke., 


fain 

S {u^B.j,—Kga^B.^i{(i)Ai)} 


and yields K^, Kg, and Kg. 

Equation (26) can bo expressed as 

Fz=- S {-KA^+KiJ^A^+Kgo>B,^-Kgo>B^-i/A^)i (3 2) 

71 jasoJl 

and yields Fg. 

In order to evaluate the remaining constants, K, K, Kg, and Fg, the coefficeint Kg is first eliminated between equations 
(23) and (24). This step is taken in order to obtam a more accurate value of the coefficient Kg, wliich is less significant under 
these conditions, and to provide a better conditioned matrix of coefficients in the computation of K, K, and Fg. Upon elimi- 
nation of the Kg term between equations (23) and (24), the following relation is obtained: 


Kg (B^B yj, -\-A^A ^) -j-Kibs {A ^5^ — B ^ A ^) — FgA.^ — 


(B ,j,B,j, -\-A tj,A,j,') — KgU)^ 




2 


(33) 
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Upon the application of the method of least squares, equation (33) becomes 


j—ia\ 

tan 

— (jsB 

J=*wt 


"AV 

tUn 

S {BfB4i-\-A^A^j{uiA^^—uB^A^} 

i—fai 


tJn 

-< 


(tin 

S {-A^}{BiB<i,-\-AfA^j 

i=«i 

J = _ 






;=ui \ 


Sri A 

' 2 \ 


which yields K 3 , K^, and F^. 


{(js^B^BaAu^AaAj, — J?sCi>*j^j ^ — oiB^A^j 

3=m\ [5r| J i 

f: ( {-Ai), I 

I=»i\ I Sri A J 


(34) 


The K 3 coefficient maj’- now be obtained from equation 
(24) w'hich, after the method of least squai‘es has been 
applied, may be expressed as: 

Ci)n 

(36) 

j“<oi 

An alternate expression for may be obtained from 
equation (23) which may be expressed as 

QK.JBff-\-KiwA^-\-K3()?B^ — «o*JB^),(wA(^)j 

(36) 

S [(w-d^)i]* 

Values of K^, which define the derivative Ct^, are important 
in the spiral mode. In order to obtain accurate values of 
Ki, aircraft motions which bring out the full effect of the 
long-period spiral mode should be analyzed to obtain 
accurate frequency-response data near <o=0. 

In the case where the solution for (7,,^ is found to involve 
ill conditioned matrices, should be eliminated between 
the real and imaginarj’' equations (25) and (26) . 

With all the K and F coefficients of the equations of motion 
determined, the lateral-stability derivatives can be computed 
from the K and F coefficients and the known aerodynamic 
parameters by use of the definitions of table I. The transfer- 
function coefficients can be calculated from the definitions 
in table II and the transfer functions and modes of lateral 
motion can be obtained from the equations in appendix A. 

SUGGESTED PROCEDURE 

In order to aid in the application of this method to the 
analysis of flight data, a suggested step-by-step procedure 
which is presented in this section has been wmrked out. An 
effort has been made to schedule the procedure so as to 
reduce the dependence of the results on the derivatives that 
can be least accurately obtained from the particular data 
being analyzed. Alternate steps are suggested where it was 
found that some particular derivative might be more 
“cen j'af“lv dotormined by onf> or th(^ other of two anproaches 


under certain conditions. No weighting of specific grou|)s 
of data is employed in the least-squares lu'oeedurc; however, 
weighting can be employed when it is corrsidered desirable 
to put more dependence on data regarded as more reliable. 

As a further demonstration of the method, two numerical 
examples have been carried out according to this in'occdure. 
In one example, presented in appendix B, the lateral-stability 
derivatives are calculated for an airplane whose transient 
responses to aileron deflection are assumed to be known. 
This example was chosen in order to gain some insight into 
the dependence of the method on the accuracy with which thi' 
frequency response can be obtained from the transient 
response. 

The other example is for the hjpotlielical rigid airplane 
whose mass and geometric parameters are listed in table 111. 
The. assumed frequency-response data of sideslip angle /3, 
roll angle <i>, yaw angle and lateral acceleration for a 
rudder-deflection input such as might be obtained from 
analysis of flight data are listed in table IV and plotted 
in figures 3 to 6. The analysis of these data has been used 
to illustrate the following suggested procedure: 

(1) Tabulate parameters and working equations 

Tabulate the airplane pai-ametei'S and any stability 

derivatives that are known or can he computed directly 
from those that are known as illustrated in table 111. 
Tabulate the least-squares equations (27) to (32) and (34; 
to (36) with any known terras on the right-hand side a.s has 
been done in equations (27a) to (32a) and (34a; to (30a) in 
table V for the example. (Since was taken to be zero 
in the example, all terms containing A'a vvere droppetl; 
however, in the example of appendix B, Cu^, was included 
and the term was determined.) 

(2) Tabulate frequency-response data 

Tabulate the amplitudes and phase angles of the four 
lateral vaiiables at the values of « to be irsed in the analysis. 
For simplicity in the example, as shown in table IV, 10 
integral values of w, evenly distributed over the I'ange of 
the data, w'ere chosen; liow'cver, in cases wlicte more accural i> 
points occur at non-integral values of w, it is advisable to 
use such points where possible. The range of values of w 
should be I'ostvicted to the rigid response. 
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TABLE II.— COEFFICIENTS OF TRANSFER FUNCTIONS 
Tlie equations for the transfer functions are given in appendix A; 
the coefficients Cn are obtained from the coefficients Cn by dividing 
by (e.g., C, =§-!)] 


C^' = l-K,Ki 

Ci' = K,+ Ki+ Kto-KiK,- KiKi- KiK,Ki 
C,' = Ki+K,Ki+ KiKio+ KiKio-KiKe-Kdu-KiKiK,- 
KJCiKi 

Cl' = KiKi+ Ki'K-- K 1 K 1 + K.KiKii- K^KiKi- KiKiK-, 

Ct,' = KiKiK^o - KiKiK-, 

Ci' = Fi(l-KiKi) 

Cl' = FiKi+FiKii-KiFi-Fi-FiKsKi-F.KiKi 

C,' = FiKiKii- F 2 K 1 + FiK2Ki-FiKi+ FiK.-F^KiKi 

Ci'=F2K2Ki-\-F2KiiK2 

C2' = F2+FiKi 

Cio' = F2K1 + F^KiKi - FiA's+ Fs Ae+ Fa Ai As+ Fj A,o 

Cii =F2 K.iK,\q~\- F F2Ki-\- FiK-i — FiAaAio“f" FiKiKi 

C,2'=Fi+KiF2 

C,a' =FaA4+ Fa Ai-FiAa As+ Fi A, + F2A9+ F, Ai Aj 

C,4' = FiKiKi - Fi A 3 A 9 + Fi AiAi+FaAi A. 

Ci5'=FaA2A7+F3A27fa 

C,i' = Ci'VFi 

C„' = Ci'VF,-VKyCi' 

Cri'=VF,C2'-VKiCi' 

C,9' = FFiCa'-FAiC7' 

C2i'=VF,CJ-VKxCi' 


TABLE III.— AIRPLANE PARAMETERS USED IN EXAMPLE 
GIVEN IN BODY OF REPORT 


(a) Known 

M=0.8 (evaluated for an altitude of 10,000 ft) 
6=22,6 ft 
5=130.0 sq ft 
p=0. 001756 slug/cu ft 
F=861.74 ft/sec 
7x = 2062 slug-ft2 
7z= 13,298 slug-ft2 
W 

r/i=— =295.03 slugs 
Ixz= 157 slug-fy 

a,=o 

(b) Computed 



= 1.5 sec 


— w 

Cl = -^=0.1121 
qb 

A6=^=0.07014 

iX 

A8=^=0.011806 
7^2=^=0.0374 per sec 


TABLE IV.— FREQUENCY-RESPONSE DATA FOR EXAMPLE 
GIVEN IN BODY OF REPORT 
(a) Sideslip, /3 


'a, 

radians 

sec 

l-l 

lirl 

deg 

\b I 

^«as — COS #j3» 

M 


1 

Sr 

sin <5/1, 

1 

0. 536 

-0.1 

0. 536179 

_ 

0. 000965 

2 

.574 

-1.1 

.674106 

- 

-.011025 

3 

.650 

-2.3 

. 649929 

- 

-. 026083 

4 

.797 

-4.1 

. 794967 

- 

-.056988 

5 

1.116 

-8.2 

1. 104368 

- 

-. 159106 

6 

2. 105 

-20.8 

1. 967417 

_ 

-. 747357 

7 

4. 388 

-109.7 

-1.462174 

— 

4. 

83763 

8 

1. 368 

-158. 1 

-1. 268821 

- 

510100 

9 

.715 

-166.1 

-. 693781 

- 

171676 

10 

.464 

-169.1 

-. 465633 

“ 

-.087720 


(b) Eoil, ^ 


63f 

radians 

sec 

111 

deg 

Idl 

~ COS$. 

ISrI 

JS4 = 

0 

Sr 

sin 0 . 


8.761 

79.0 

1.676642 


8. 598903 


4. 624 

67.7 

1. 754686 


4. 278325 


3. 432 

57.4 

1.850156 


2. 891147 


3. Ill 

47.6 

2.096761 


2. 297575 


3. 459 

37.1 

2. 760823 


a 084203 

6 

5.446 

18.9 

5. 153317 


1. 762687 


9.719 

-74.6 

2. 578202 

— 

9.370572 


2.729 

-126.9 

-1. 638908 

— 

2. 181509 


1.298 

-138.2 

-. 966715 


-. 865863 

1 “ 

.780 

-143.9 

-. 630341 


-.459096 


(0) Yaw, li' 


radians 

sec 


t 

~Sr 


*4, . 
deg 

i^l 

Ir 1 cos^^ 
1 1 



0. 223 

163.4 

-0. 213722 

0.063714 



.494 

174.4 

-. 491634 

.048463 



.613 

174.4 

-.610140 

. 059934 



.773 

172.8 

-.767184 

.096776 


1.096 

168.9 

-1. 076250 

. 210567 

6 

2.079 

156.2 

-X. 903129 

.837405 

7 


67,3 

1. 656940 

3. 966906 

8 

1. 358 

18.7 

1. 285897 

. 436504 

9 



10.7 

. 698169 

. 13J797 

10 


.462 

7.1 

. 458498 

.057029 


(d) Lateral acceleration, a. 


radians 

sec 


^fl 

ft/sec» 



deg 


radian 


1 

2 

3 

4 

5 


107. 67 
121. 70 
149. 83 
204. 00 
322.10 


179.7 

178.2 

176.5 
174.0 

169.6 


dy 


cos 


ft/scc^ 

radian 


-107.6736 
-121. 6356 
-149. 6411 
-202 8943 
-316. 7986 



Ory 

It 


sin <f>a 


ft/sepg 

radian 


0. 564215 
3. 928440 
9. 277399 
21. 205800 
58. 143795 


6 

7 

8 
9 

10 


691. 42 
1628. 37 
587. 34 
350. 66 
254. 54 


156.5 

67.3 
18.6 

10.4 
7.2 


-634.3160 
627. 9656 
556. 6625 
344. 9205 
252. 5535 


275. 1^748 
1502. 418071 
187. 337711 
63. 200245 
31. 756598 



Figure 3. — Frequency response in sideslip due to rudder-deflection in- 
put for rigid, high-speed airplane of example given in body of report. 
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Figure 4. — Frequency response in roll due to rudder-deflection input 
for rigid, high-speed airplane of example given in body of report. 





Frequency of rudder-deflection input, rodians/sec 

Figure 5. — Frequency response in yaw due to rudder-defleclion input 
for rigid, high-speed airplane of example given in body of report. 

TABLE V.— MODIFIED FORM OF EQUATIONS USED IN EXAMPLE GIVEN IN BODY OF REPORT 




(27a) I 
(28a) 
(29a) 



Urn / 

nF^=-j: ( 

j=m \ 

K,A,+^^ 

J 

(30a) 


f) (o^am-b,),' 

— \ 

'2 i.i»^B^-K,b>^B;j,{-B^)j 

i 

Wr 

2j i~S0)j(blA^)j 
-/=«! 

/— <dl 

2 Ko>A,hr 

i=«i _ 

T 

11 

o 

^®*wi 

a>B 

(w’^Bi^ — Kib>~Bf)j(uA^)j 
tf=wi ^ 

^ (31a) 1 




{B^B^-\-A^^){caA^B^ — (cB^A,f)f 


i—m 

■ «w Ww 

2 {<^A^B^—(/>B^A^/^B^B^AA^‘i^)j 2] {—A^}{B^^-\-Afjx\^} 

j=m j=u>i 


(32a) : 


S (.B0B^+A0A^)j(—A^)j 
/=«! 


y=«i 

^ ) {<t)A(f,B^ — <aB,^A^}^ — 

\ \8t\ } j 

i^rl } s 

j="l \ \8r\ /j 

tOn 

S iB’2+u’‘A^—Ksa’‘A^+KiuB,^—KiAfi) }iuB^)j 


-m—m 

( A.}^j ((liA^B^ (xi Bt*iA\^j 


2 

f = Wl 


^ K, 
IF,J 




(:Ma) 






S \W,Y 

S (-^3-Bs -{-KibiA ^ -f-if s«w®5 ^ — b?B ^j{uA^] 


)=a, 


(35a) 


(36a) 


Yaw phase angle, 
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PiGUHE 6. — Frequency response in lateral acceleration due to rudder- 

deflection input for rigid, high-speed airplane of example given in 

body of report. 

(3) Compute vector components of data 

Compute the vector components from the data listed in 
step (2) and tabulate them as shown in table IV. It is 
advisable to use at least 5-place tables of sine and cosine 
values in this computation. Additional columns of 
0)5,3, . . . and <AA^, . . . vtU be 

useful in subsequent operations but are not listed for the 
example. 

(4) Compute Ki 

Compute iifi from equation (27a). (See table VI (a).) 
In this step and in subsequent steps, summation of the 
corresponding elements of the indicated products of the data 
columns evaluates the matrices of coefficients and knowns. 
The resulting simultaneous equations are readily solved by 
any standard method. 

(4a) Alternate step 4 

Compute Ki from equation (29a). (See table VI (o).) 

(5) Compute 5) 

Compute Fi from equation (28a) with the values of Ki 
from step (4) . (See table VI (b) .) 

(5a) Alternate step 5 

Compute from equation (30a) with the values of Ki 
from step (4a) . (See table VI (d) .) 

(6) Compute K, and Kio 

Compute Kj and Kw from equation (31a). (See table 
VI (e).) 

(7) Compute Fz 

Compute Fz from equation (32a) with the values of Ki 
and Kio from step (6) . (See table VI (f) .) 

(8) Compute Kz, Ki, and Fz 

Compute Kz, Ki, and Fz from equation (34a) . Because 
of the combination of two equations to obtain equation (34) 
in the method, the combination of data columns is more 
complicated in this step than in others. (See table VI (g)..) 

(9) Compute Kz 

Compute Kz by one of the following methods: 

(a) If Fz has not been eliminated, calculate Kz from 
equation (35a) with the values of Kz, Ki, and Fz from step 
(8). (See table VI (h).) 

(b) If Fz has been eliminated, calculate Kz from equa- 
tion (36a). (See table VI (i) .) 

(10) Compute stability derivatives 

Compute the lateral-stability derivatives from the 
relationships in table I and the now complete set of values 
of K and F. (The computed and known values of the K and 
F coefficients and the stability derivatives for the example 
are given in tables VII and VIII, respectively.) 


(11) Compute transfer-function coefficients 

Compute the transfer-function coefficients from the 
relationships given in table II and the values of K and F. 
(The computed and known values of the transfer-function 
coefficients for the example are given in table IX.) 

(12) Compute the modes 

The modes of lateral motion may now be computed, if 
desired, by use of the equations in appendix A. (The values 
of the lateral modes for the example are given in table X). 

(13) Accm-acy checli 

When the true derivatives are unknown, the check on the 
computations is obtained in the follovdng manner: Calculate 
the frequency response from the equations in appendix A 
and the computed coefficients of step (11). Compare tliese 
values with the data of step (2). 

DISCUSSION 

DATA REQUIREMENTS 

The method of this report is based on the assumption that 
equations (1), (2), and (3) represent the lateral motions of 
the airplane. These linear differential equations are satisfied 
by small variations in roll, yaw, and sideslip but will not 
necessarily hold for large changes in these variables or in 
regions of flight where nonlinearities occur. It is therefore 
desirable to obtain data for transfonnation to the frequency 
plane from small motions recorded within ranges of Mach 
number and other variables in which the coefficients are 
constant; also, the range of frequency-response data to be 
used must be limited to values for which the aircraft acts as a 
rigid body and higher frequency instrument inaccui’acies are 
avoided. Thus, the assumptions of the method will be most 
nearly in accord with the true physical nature of the problem. 
In any case, experience has showm that the derivatives 
obtained by this method of analysis will accurately reproduce 
the frequency response from which they were extracted. 

Since the present method requires flight data in frequency- 
response form, either the forced-oscillation technique or the 
technique of oscillating the rudder or aileron at various 
frequencies and measuring the steady-state response in /3, 
D^, D\j/, and must be employed or the transformation of 
a tr'nsient response from the time plane to the frequency 
plane must be made by an appropriate method. A compre- 
hensive survey of methods for effecting this transformation 
from time to frequency plane is given in reference 8 which 
evaluates the Fourier integral method as w'ell as the use of 
devices such as the Fourier synthesizer, Coradi harmonic 
analyzer, and IBM machines for obtaining frequency re- 
sponse from transient-response data. 

If, as is usually the case, the transient responses in roll 
and yaw are measured in the form of time histories of 
and D\l/ rather than of 4> and 4', it is recommended that these 
time histories of Z>0 and Df be converted to frequency 
response initially and then the frequency responses of (j) and 
4f be computed from the frequency responses of D<j> and D\]/ 
by use of the relations 


D(i> 


i. 


D4i 


t 

Ir 

= 03 

5. 



= 03 

Sr 


and 
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TABLE VI.— LEAST-SQUARES SOLUTION OF EQUATIONS 
FOR DATA OF EXAMPLE GIVEN IN BODY OF REPORT 

I i 

' (a) Ki, equation (27a) 


S (-£oA^-uA#-l-X2Bi,)S*=7.378 
10 

17.562 

W=1 

17.562A’i=7.378 
Xi= 0.420 


(b) Fi, equation (28a) 


10 

^ ) (XiAjff — — wjB(j — a)B^) = 1.061 
(0= L 

10Fi= 1.061 
Fi= 0.1061 


(c) Ki from lateral-acceleration data, equation (29a) 


10 

S (-US^)2= 13,042,324 
10 

S (- VBf) = 5,567,902 

U=1 

13,042,000A'i=5,568,000 
Ai= 0.427 


(d) Fi from lateral-aocelefation data, equation (30a) 


Z:(aV1.+^“)=I.03S069 

1 OF, = 1.035 
F, = 0.1035 


(e) A 7 and ATio, equation (31a) 


2 (-5^)2=17.562397 

Ci>asl 

10 

s (-D/j)«A# = 44.47131 

(d^l 

10 

S (wA#)2= 473.928987 

10 

s A 8 co=B*) (- Bfl) =856.325838 


TABLE VL— LEAST-SQUARES SOLUTION OF EQUATIONS 
FOR DATA OF EXAMPLE GIVEN IN BODY OP REPORT— 
Continued 


(g) K 3 , Ki, and Fj equation (34a) 


S (5^5^+ -4tf.Q)2= 371.762380 

(Kl^l 

10 

A,i-1^) (wA^B^ — 0)5*.!^) = — 3617.077695 

til=l 

10 

2D + 23,184790 

W*1 

10 

XJ (uA*B,f -uBjA;.)2=36234.43887.S 

10 

S (uA*B,t-ioB,Af)f- A#)= -209.219323 
01 = 1 

10 

^ (-A^)2= 1 1.122853 

W=1 


+ )(B^B; + AsA#j= 33187.734509 | 

u=l\ 1*1 / 

10 / |T 1 J\ I 

[ w^B^B^-pw^A^A,^ — Aati,2j=-j jicoA^B^ — taB^A^)= — 322344.219155 

u=l \ 1*1 / 

u 2 B#B,f -p 0,2 A, A,, - AW i V _ J = 2422.633922 

m=i\ 1*1 / 

371.76238A3-3617.67769A',-P23.18479Fj = 331S7.73451 
- 3617.67769X3+ 35234.43888 A'l - 209.21 932F3 = - 322344.2 1 91 5 
23. 18479 A'j - 2O9.21032A\+ 1 1 .1 2286Fj= 2422.633922 
X3=138.272 7G=5.212 >2=27.636 


(h) Ee, equation (35a) 


S («D#)2=811.55027 
10 

^ / (wB^) KibuD^ — = 235.1 17»592 

W = 1 

813.Giv6=235.1 

7^4=0.290 

(i) Alternate means of obtaining Jv®, equation (3Ga) 

10 

2 C«A^)2 = 473,9280S3 

W = 1 

10 

2 = 137.774050 

(ii=i 

473.97vj= 137.77 
A'6=0.291 


2 («“B#-Xg<d2B*)(£oA;A)=23S6.79I416 

<d=l 

As= 0.011806 7f:9=0 

17.567G+ 44.47Xio= 856.3 
44.477i7-P 473 . 93 X 10 = 2356.8 
X,= 47.443 
X,o= 0.5217 


(f) Fs, equation (32a) 


10Fs= (— XyAfl-pZiLs^^A^— Xio£s/J,f—cAA^) = — 251.87 

co = l 


Fs= -25.187 


TABLE VII.— COMPARISON OF COMPUTED AND KNOWN 
VALUES OF COEFFICIENTS OF EQUATIONS OF LATERAL 
MOTION FOR EXAMPLE GIVEN IN BODY OF REPORT 


CucOicicnt, 

Computed from 1 
data . 1 

Known 

value 

Ki 

f *0.42054 
1 to. 4200 1 

[ 0.427 

IC, 

138.272145 ' 

138. 245 

K, 

5. 2125GS 

5. 21 

Ki 

0.290 

0.3017 

K-, 

47. 438G2.5 j 

47. 41 

Kii 

0*521457 1 

0. 5272 

F, 

/ *0.10009 

1 ‘0. 1035 1 

[ 0. 104 

Ft 

27.030187 

27.05 

Fi 1 

-25. 187091 

-25. 22 


•From sideslip data. 
i"From lateral-aceeleradon data. 

















DETERMINATION OF LATERAL STABILITY DERIVATIVES AND COEFFICIENTS FROM FREQUENCY-RESPONSE DATA 11 


TABLE VIII.— COMPARISON OF COMPUTED AND KNOWN 
STABILITY DERIVATIVES FOR EXAMPLE GIVEN IN 
BODY OF REPORT 


Stability 

derivative 

Computed 

Known 

Cyfi 

-1.259 

-1.28 

Ci, 

-0. 1490 

-0. 149 

C'jp 

-0. 4280 

-0.428 

Cl, 

0. 028S9 

0.0248 

GnP 

0. 3292 

0. 329 


-0. 2732 

-0, 279 

Cl-, 

0.3183 

0.312 

c,,’ 

0. 02979 

‘ 0. 0298 

r 

-0. 1748 

-0. 17o 


TABLE IX.— COMPARISON OF COMPUTED AND KNOWN 
COEFFICIENTS OF TRANSFER FUNCTIONS FOR EXAMPLE 
GIVEN IN BODY OF REPORT 


Coefficient 

Computfiti from 


data (least-squaros 
value; 

Known 

Co' 

0.999101 

0.999101 

Cl' 

0. 150763 

6. 160254 

Cl' 

50. 934260 

50. 972855 

Cl' 

253. 483361 

253. 2143 

Ct' 

2. 193801 

'2.190861 

C(' 

0. 105995 

0. 103907 


25. 488778 

25. 48986 

Cl' 

132.644801 

132. 64415 

Ci' 

C,' 

0.267223 

0. 200611 

25. 718442 

25, 729749 

Cm' 

3.636361 

3,962748 

Cii' 

-2174. 851610 

-2178. 769199 

C)/ 

Cl3' 

-24. 860818 

-24. 893564 

-136.884812 

-137. 264852 

On' 

-28.978823 

-30. 417643 

Cii' 

-81. 654072 

-81. 36940 

89. 111766 

89. 62096 

Cit' 

C,,' 

609. 187633 

514. 350717 

-4822, 026985 

-4815.419162 

Cjg' 

Cso' 

-26127. 620554 

-26138. 324577 

97. 287337 

LOO. 542037 


TABLE X.— COMPUTED LATERAL MODES FOR EXAMPLE 
GIVEN IN BODY OF REPORT 

[CharaotorisU'o equation is equation (A16) of appendi.t A] 

Damping-in-roll mode, Ai — 5 . 393 

Spiral mode, Xj - 0 . 008668 

Dutch roll mode, —i±in- - —0. 381±«(6. 84) 

Flight data are usually measured about the body axes; 
before they can be used, however, cither they must be con- 
verted to a form corresponding to the stability axes or the 
development must be restated in terms of body axes. The 
following equations are presented to be used in transforming 
to the stability axes: 

4>=^i, cos a-f-^6 sin a 

^=^6 cos a— 06 sin a 

where a is the angle of attack (angle between the body A'’-a.\is 
and the stability JY-axis) and the subscript b refers to the 
body axes. In addition, the sideslip-vane measurements 
can be converted to the corresponding values at the center 
of gravity by means of the formula 



where the subscript v refers to the vane and is the distance 
from the center of gravit}' to the sideslip vane measured 
along the body A"-axis. 

If the lateral acceleration a, as recorded by an accelerom- 
eter is used, it must be collected for any components of 
angular velocity and acceleration present in such records 
because of the displacement of the accelerometer from the 
airplane center of gravity. 


SENSITIVITY OP THE METHOD AND POSSIBLE SIMPLIFICATIONS 


Even though the foregoing general rules are observed, the 
question arises as to how the derivatives obtained compare 
with the actual derivatives of the airplane and how sensitive 
these derivatives are to inaccuracies in the frequency- 
response data. These inaccuracies stem from several sources 
such as measurement of the transient motions, transformation 
from time to frequency plane, and reading errors. Care 
should be exercised to keep these inaccuracies to a minimum. 

No attempt is made to present a comprehensive error 
analysis in this report, but several general remarks appear 
to be in order. A comparison of the values given in tables 
VII, VIII, and IX for the coefficients of the equations of 
lateral motion, the lateral-stability derivatives, and the 
transfer-function coefficients indicates essentiallj’’ no errors 
due to the method itself. It has been found, however, that 
the derived values of Cy^ and Cy^^ are sensitive to slight 
variations in the quantities A^, Af, and (see eqs. (27) 
and (28)). This sensitivity is due to the fact that the air- 
plane angle of yaw is very nearly equal to the negative of 
the sideslip angle and hence accurate determination of the 
lateral acceleration from measurements of 0 and /3 is diffi- 
cult. This difficulty can be obviated by a direct measure- 
ment of the lateral acceleration with an accelerometer and 
more accurate values determined for Cy^ and C'va by means 
of equations (29) and (30). 

As brought out in appendix B, the accuracy with which 
the values of the derivatives Ci^, C„^, and Cy^ can be deter- 


mined from frequency-response data appears to be depend- 
ent upon the inclusion of data for frequencies near the natural 
frequency of the airplane. This fact is brought out further 
by some unpublished results that indicate that derivatives 
such as Cn^ and Cy^ are very sensitive to variations of the 
frequency responses near the airplane natural frequency. 
In the case of On^ the conclusion is obvious smee this deriva- 
tive is laiown to contribute largely to the airplane Dutch 
roll damping and it is also loiown that the peak value of the 


frequency-response amplitude ratio, say 


15 


, is highly de- 


pendent upon the airplane Dutch roll damping. It appears 
thei-efore that an accurate representation of the frequency 
responses at frequencies near the airplane natural frequency 
(for instance ±50 percent) is probably very important in 
the determination of the lateral-stability derivatives. 

In general, when the number of unknown derivatives can 
be reduced, the remaining derivatives can be determined 
with greater reliability. For this reason,, several possible 
simplifications are presented for the calculation of the 
derivatives (7;^, and Ci^^. The airplane natural 
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frequency is, for small values of the product of inertia, well 
approximated by the expression 


Wjv‘ 


'S 


relation. The value of 


to zero frequency 


erivative Cn^ subse- 


at ca=0 obtained by this extra- 


If the natural frequency is apparent from the available 
transients and the inertia parameter is Imown, a vei-y good 
value of Cnp can generallj’- be obtained from this expression. 

I a 
Or 

may afford a m.eans of obtaining the c 
quent to determ.ination of from the preceding frequencj- 

i 

Sr 

polation is not indicative of the actual static sensitivity, 
but rather of the apparent static sensitivity which includes 
the veiT lightly dam.ped spiral mode. This apparent steadv- 

state response is approximately equal to the ratio 

from which can be estimated if (7,,^ is known. 

The actual static values of some of the lateral variables 
may be found theoretically by taking the ratio of the integrals 
of the time histories of the output and input; for exajnple. 


M 0 


n CO 

.J> 

" fco 

J. 


Probably in certain cases these integrals could be evaluated 
accurate^ enough to yield a good approximation to the actual 
static sensitivities. Some of the stability derivatives 
influence greatly the values of the transfer functions for 
u=0, and use of accurate zero frequency data would 
undoubtedlj^ lead to better estimates of these stability 
derivatives. 

The transfer function which relates the rolling velocity to 
aileron deflection is well approximated (particularly at low 
frequencies) by the expression 





(?2 -I tan-1 


to «=0, the parameter and hence 


, • '<l> 

amplitude ratio 

can be obtained once Ki has been determined. x\s was the 

|/3| 


case for 


Sr 


, this extrapolation does not yield thi' true value* 


of 


at £ 0 = 0 , which is actually zero, but gives the apparc'iit 


value due to the lightly damped spiral mode. 


CONCLUDING REMARKS 

A method has been presented for extracting lateral-stability 
derivatives from frequem-y-response data which have 
been derived from aircraft transient responses to arbitrar,\- 
control inputs. In order to demonstrate the use of the 
method, the lateral-stability derivatives have been calcu- 
lated for two hypothetical aii'planes for which the frequency- 
response data or transient-response flight data were assumcil 
to be known. Simplifications were proposed for obtaining 
certain of the derivatives, and although no error analysis is 
presented in this report, some general observations can l>e 
made concerning the sensitivity of the method to inaccuracies 
in the original data. 

In view of the limited experience in the determination of 
lateral-stability derivatives from flight data, and particulai’ly 
in using the present method, it appeam that an inve.stigation 
should be made to determine the effects of some of the orrons 
inherent in analyses of transient responses, such as recording 
accuracy, reading accuracy, and so forth, on the dei’ivetl 
frequency responses and hence on the computed latei*al- 
stability derivatives. It is noted here that, although the 
methods presented hi this report foi* the extraction of lateral- 
stability derivatives from frequency-response data are 
theoretically correct and matheraatiaxlly sound, the accuracy 
with which the lateral-stability derivatives of an airplane 
can be determined by this method is directly depe-iuh'nt upon 
the accuracy with which the frequency responses of the air- 
plane are Icnown. Also, some consideration should be given 
to the types of inputs likely to afford good frequency-i’t*- 
sponse data and to testhig procedures whereby some of the 
stability derivatives that are most difficult to determine may 
be accurately obtained. In addition, accurate frequency- 
response data for specific ranges of frequency may yield 
more accurately certain of the stability derivatives and this 
possibility should be mv^estigated. 


Therefore, (,w) = — tan“‘ ^ from which Kt, and subse- 

quently can be determined. Also, by extrapolating the 


Langley Aeronautical Laboratory, 

National Advisory Committee for AEito.vAUTrcs, 
Langley Field, Va., February 16, 1954- 



APPENDIX A 


DEVELOPMENT OF EQUATIONS FOR TRANSFER FUNCTIONS, MODES, AND FREQUENCY RESPONSE FOR LATERAL MOTION 


In this appendix the transfer functions are derived from the 
equations of lateral motion and the modes and frequency 
response are expressed in terms of the transfer-function 
coefficients. 

EQUATIONS OF MOTION 

The motions of a rigid aircraft resulting from a rudder- 
deflection input are assumed to be expressed by the following 
thiee standard lineaiized equations of lateral motion: 

(Al) 

ifa/3+ (Kj:>^+KM=FMt) (A2) 

{KJD^+K^) <t> + (D^+K,,B) >p=FMt) (AS) 


where the K and F coefficients are defined in table I. The 
axes and the sign conventions employed are shown in figure 1. 
On applying the Laplace transformation 




‘dt 


(A4) 


(where a=o-ff«) to both sides of equations (Al), (A2), and 
(A3) and assuming the initial conditions to be zero, there is 
obtained 


— Ki s 

s'^+K,s -K,s^-K,s 

l-K, -K^^-K,s s^+K, 

TRANSFER FUNCTIONS 



'13(g)' 


'W)' 


<^>(s) 


F^SrCs) - 




.BsSr(s)^ 


) 


(A& 


(A7). 


By solving the three simultaneous equations (A.5), (A6), 
and (A7) for p{s) there is obtained 

(C,'s^+C,'s^+ C'a'sH Ci's)m 

= (C-/s^+C,'s^+C/s^+Cs's)8M (A8) 

where the On coefficients are defined in table II. Dividing 
equation (A8) through by Co' and using the unprimed 
coefficients C„ to designate the ratio of the primed coefficient 

to Co'( e. g., Ci=^ and 02 =-^, j results in the following 

transfer function for /3 due to a rudder-deflection input 5r: 

/3 Cos'+Cos'H-as-l-Cs 


hr S'*-l-CiS®-l-C2S^-t-C3S-l-C4 


(A9) 


Determined in a similar manner, the transfer function for 
^ is 

4> CpS^-f-Ciog-l-Cii 


Sr S^-|-CiS®-f‘C 2 S^-l-f%S-|-C 4 
and the transfer function for ^ is 

^ Ci8S^-hCi3S^-|-Ci4g-|-Ci6 


5, g(g^H-Clg®-t-C2gHC3g-fC4) 


(AlO) 


(All) 


The lateral-acceleration transfer function may be developed 
from the sideslip transfer function in the following manner: 
The lateral acceleration is usually defined as 



ay=V(D^+Dyl/—K2<t>) 

(A12) 

Substituting equation (Al) into (Al2) yields 



ay=VFMt)-VK,^ 

(A13) 

Therefore, 

^=VF,~I^ VK, 

Sr Sr{S) 

(A14) 

or 

ay C'iBg^-t-C'i7g®-f-6'igg®-|-6''i9g-l-C20 

(A15) 


Sr S^+O^S^+C^S^ + C^S + C, 


MODES 

The characteristic equation of lateral motion is obtained 
by expanding the determinant of the coefficients of equations 
(Al), (A2), and (A3) and is 

g(g*-fC,g'-t-C2g'+C3g+C4) = 0 (A16) 

This equation can be factored as follows : 

s(g+ Xi) (g-|- Xs) (g+ ^ — 11?) (s +?+'!»/) = 0 (Al 7) 

where the root g=0 indicates that the aucraft is insensitive 
to azimuth, the root Xi is the damping-in-roll mode, the root 
X2 is the spiral mode, and the root — is the Dutch i-oll 
mode, an oscillatory motion composed of roll, yaw, and 
sideslip. The modes of this airplane are indicated in table X. 
An approximate expression for X2 which. gives excellent re- 
sults is 

X3=&=S (A18) 

^3 ^3 


FREQUENCY RESPONSE 

The frequency response which is the steady-state response 
to a sinusoidal input consists of an amplitude-ratio and phase- 
angle relation between input and response for various values 
of «. It is determined by substituting s=iu into equations 
(A9), (AlO), (All), and (Al5), respectively, for |3, <j>, yp, and 
ay and is a complex expression. The amplitude ratio is the 
square root of the sum of the squares of the real and imagi- 
nary parts of the complex expression and the phase angle 
$ is the ai’c tangent of the ratio of the imaginarj’- part to the 
real part. 

The amplitude ratio for sideslip angle ^ is 

IL I (..a) 

Sr\ V(«^-t?2Co‘*+C4)H(C'3W-C7lCoy ^ 

and the phase-angle relation between ^ and 5, is 


$,35 =tan“* 


/ C7C0 — C5W® 



0 ) 


CsU — CiW® \ 


(A20) 
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1* 


(A21) 


The amplitude ratio for roll ^ is 

■^1= / (f^ii~^o t^ y+(C''ioto) 

8rl \ CjcoH C,f + {CiW-i 

and the phase-angle relation between <j) and 5, is 
Tbe amplitude ratio for yaw xj/ is 

^[_ / (<^15~^13t«>^)^+(C'i4fa) — CiaCO^)^ /■AOQ\ 

5,1 \ 


and the phase-angle relation between x{/ and 6^ is 

(A34J 

Tlie amplitude ratio for lateral acceleration is 

(A25j 


Sr 




'■ — f^l8a>^d~(^2o)^4~(C'i9tO— C'i7td*)* 


and the phase-angle relation between a,j ami 6, is 

=tan- ‘ (■■ )-t an- ■ ( 

Vt-ifico* — Ci8^ ~r^- 20/ V&j^—Caw^+Cv 


(A2G^ 



APPEPTOIX B 


DETERMINATION OP LATERAL-STABILITY DERIVATIVES BY USE OF FREQUENCY RESPONSES DERIVED PROM CALCULATED 

TRANSIENT RESPONSES TO AILERON DEFLECTION 


Transient responses in roll, yaw, and sideslip to a square- 
pulse aileron deflection were calculated for the airplane 
whose mass and aerodynamic characteristics are given in 
table XI. These transients and assumed aileron time his- 
tories are presented in figure 7. 

TABLE XI.— PARAMETERS OF AIRPLANE USED AS 
EXAMPLE IN APPENDIX B 


Altitude, ft f 30, 000 

Wing loading, Ib/sq ft 65 

V, ft/sec 797 

6, ft 28 

Cl - 0.23 

HI) 80. 7 

Kx^ - 0.0097 

AV - - 0.0513 

Kxz 0.00145 

C(p per radian —0.40 

per radian 0. 08 

per radian - —0.02 

Cn, per radian —0.40 

per radian —1.0 

Cn^ per radian 0. 25 

Ci 0 per radian — 0. 126 

Cij per radian —0.11 

"a 



Time, sec 

(a) Rolling velocity. 

Figure 7. — Calculation of transient responses to a square-pulse aileron 
deflection for the airplane defined in table XI. 

The frequency responses D<l>l8a, Dfjha, and were ob- 
tained from analysis of the calculated transient responses 
D<j>, Vyj/, and /3 to the assumed aileron deflection by use of 
the equation 

e-“Fif)dt 

h 

e~^‘q(f)dt 



Figure 7. — Continued. 


32x10''^ 



(c) Angle of sideslip. 
Figure 7. — Concluded. 


curves and thus perfoi-ming analytically the required inte- 
grations from 0 to <» in finite time intervals as described in 
reference 8. These transforms were then evaluated for a 
number of values of iu (where s=ico) and frequency-response 
data were obtained. These data have been plotted in 
figure 8 where they are represented by the triangular test 
points. 

For purposes of comparison, the frequency responses for 
the airplane were calculated from the basic rrrass and aero- 
dynamic characteristics and are also presented in figure 8 in 
the form of curves. The agreement between these two sets 
of data is seen to be e-Kcellent. 

In addition, the frequency responses were obtained from 
the calculated transients by use of IBM equipment. The 
Fourier transform of each transient was obtained from a 
numerical evaluation, by means of the IBM equipment, of 
the integral 


The Laplace transforms of the responses were obtained by 
fitting various functions, such as polynomials or trigono- 
metric functions, to finite sections of the respective output 


fe-^-*F{t)dt 
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Figtibe 8 . — Frequency responses due to aileron-deflectioii input for 
the airplane described in table XI, 




w, rodions/sec 


(b) Pliase angle 
Figttee 8. — Continued. 


Figure 8. — Continued. 



Figure 8. — Continued. 


between the limits of t=0 and f=6.9 seconds to which was 
added an end correction based on an analytical expression 
of the transients from i=6.9 seconds to infinity. A time 
interval At =0.05 second was used in the numerical integra- 
tions. The frequency responses thus obtained are shown in 
figure 8 as square points and appear to be in excellent agree- 
ment with the frequency responses calculated from the basic 
mass and aerodynamic characteristics and with those derived 
from the transients by the curve-fitting method. 

By use of equations (27), (31), (33), (35), and (36) of the 
body of the report and the frequency-response data derived 
from the transients by both methods, the coefiicients Ki, Kz, 
Ki, Kz, <? 2 , Ki, Kz, and isTio were obtained, and hence the 


stability derivatives C,.^, Ci^, Ci^, Ci^, C„^, and 

These results are presented in table XII along with 
the values assumed for the calculation of the transii'iit 
motions. Results are also presented for the case whore' 
equations (1), (2), and (3) of the body of the report wore 
not separated into real and imaginary parts but were troati'd 
as vector equations and the method of least square.s a|>i))iod 
in a manner similar to that described in reference 4. TIu* 
amount of work required in the least-s<iuares process as 
shown in appendix C is less when these equations are treated 
as vector equations, but a limitation is introduced in that 
a parameter such as K wliioh is difficult to extract with 
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Figure 8. — Continued. 
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Figure 8. — Concluded. 
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good accuracy cannot be eliminated as it was in tbe body 
of the report when the equations were separated. 

A comparison of the values of the derivatives based on 
frequency responses obtained by the curve-fitting method 
and obtained from the two methods of treating the equa- 
tions of motion indicates very good agreement for aU the 
parameters except Ot^. The poor agreement for Ci^ is felt 


TABLE XII.— COMPARISON BETWEEN STABILITY DERIV- 
ATIVES USED IN CALCULATION OF TRANSIENTS AND 
THOSE DERIVED FROM TRANSIENTS FOR EXAMPLE IN 
APPENDIX B 


Derivative 

Assumed 

value 

Transient analysis, curve- 
fitting method 

Transient analysis, IBM 
equipment 

From eqs. of 
present 
report 

From eqs. of 
motion as 
vector eqs. 

From eqs. of 
present 
report 

From eqs. of 
motion as 
vector eqs. 

■ ^‘v 

-0.400 

-0. 378 

-0.397 

-0.388 

-0. 391 


0.080 

0.0882 

0 . no 

0. 130 

0. 099 


-0. 126 

-0. 127 

-0. 134 

-0, 129 

-0.128 

< 

-0. 100 

-0.100 

-0. 100 

-0. 098 

-0.099 


-0.020 

-0.020 

-0.021 

-0.029 

-0. 025 

c/ 

-0.400 

-0. 382 

-0. 397 

-0. 786 

-0. 646 

Cn- 

. 0.250 

0. 253 

0.257 

0.213 

0.258 

9', 

-1/000 

-1. 050 

-1. 003 

-0.600 

1. 260 


not to be significant since Ci^ generally makes a negligible 
contribution to the lateral motions of an airplane and is 
important only in the lightly damped spiral mode of motion 
which is very difficult to define accurately from a transient 
response. 

A similar comparison of the two sets of values based on 
frequency responses obtained with the IBM equipment 
shows good agreement except for the derivatives Ci^, On., 
and The poor agreement between the initially as- 

sumed values of (7?^, and Cy^ and those calculated from 
the frequency responses obtained with the IBM equip- 
ment is probably due partly to the fact that data were not 
obtained for frequencies near the natural frequency of the 
airplane with the particular routine used for the IBM 
equipment and partly to small inaccuracies that exist in 
the data at the frequencies that were included. These 
small inaccuracies undoubtedly could have been reduced 
somewhat by the use of a smaller time interval in the numer- 
ical integrations performed by the IBM equipment. 

The use of the frequency-response data which were ob- 
tained from the transient responses by the curve-fitting 
method and which included frequencies near the natural 
frequency of the airplane resulted in good agreement for 
these derivatives. Therefore, it appears that further investi- 
gation into the effects of the choice of the frequency 
range to be used in extracting stability derivatives from the 
frequency-response data is warranted. 





APPENDIX C 

LEAST -SQUARES VECTOR SOLUTION OF THE EQUATIONS OF LATERAL MOTION 


Some of the results shown in table XII were obtained bj^ 
applying the method of reference 4 to the lateral equations 
of motion (appendix B). The important details of this 
method are illustrated in this appendix by appljdng them to 
equation (1) which is 


[Sr(0] 


(Cl) 


If the flight data are expressed in frequency-response form, 
that is for every frequency wj there is an amplitude ratio and 
a phase-angle relationship, then the lateral variables can be 
considered to be constant-amplitude vectors rotating in the 
complex plane with fixed phase relationships for any given 
freqiiency. Equation (Cl) can be expressed as 


O'-"..".. •■•“j 

(C2) 

where airplane frequency re- 

sponses expressed in vector form, and the residual Vj is a 
vector at each frequency. In this treatment K> is assumed 
to be known and the most probable values of the unknowns 
Ki and Fy are such as to make the sum of the squares of the 

residual a minimum; that is, {ojY is equal to a minimum. 

The squai'e of the residual vector is obtained from the scalai 
or dot product of the two vectors, nanieljq 




(C3) 


The values of Ky and Fy which satisfy the condition that S be a minimum are obtained from a simultaneous 

J—bn 

solution of the equations 


5jr2fe.»,)=0 

These equations after expansion are written in matrix form a.s follows: 
(a) For the sideslipping motion 


j—<»n 


(b) For the rolling motion 


S (~Afi)j 


/I fll\ ^ 

^(m 

J=wi \|0r|/j _ 


Fy 

Ky 




^ — (aB ^ — 

U>n 

'y~\ -j- BgB,^\ j 

j=m 


n 

2 


i=«i 

(an 


ZC-A3)i 

I 

j*«l 



(an 

2 uj(~A^B^-\-A^fi)j 

j=<ai 


can 

2 

i=ui 

2 c^i(ApB^—A^Bff)j 

J=(ai 


2 

}=m 




21] — 


2 (— w5lf)j 

3—a\ 

tan 


]2 — a>j\A^^-\-B^B^j 


Q>i> / \th 

) 

I Sr A 




K, 

lU 

lu 


>■ = 


2 w/(— /!C5.fi^)/ 




(C4) 

(05) 


2j &>/(AflA0+5(3B, — K-^fiA^—K^jiB^)j 


j=tn 


2 K^(>)j\A^B^—B^A^) 


tan 

z 

j«u>l 
ta» 


S (ii/(A^B^—A^B^)] 

i=vit 


t(.'7) 


3S 
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(c) For the yawing motion 


n 




Z 

Un 

Z (,‘>>B^)j 


F, 





j—ui 

j—Ul 





z(|[) 

Un 

z 

Uj( — A^^-]-A^B0)j 

(>>j{A^Bj,—A^B^} 


K, 

j—m 


J=m \10r| / j 




H 

V 


0)n 

So 

oj(—AffB^+A^Bff)j 


,!;(“■ ID 

Z — (>ij\A^^-\-B,^B^j 


K, 

j=ui 




J-ai \ /j 

§—m 



6in 

■ j—on 

(t>fi 

Z! 

o)j(AfiB^—A^B^j 

<an 

z 

j=ui 

—uJ‘(A^A^-\-B,jyB4f)j 

\6r\/3 


A'lo 


g>ft 

j=o>i 


tii>n 

i*wi 

«iL 

t*>/( — A^By^-\‘A^B,^j 

J—m 

Uh 

KiUj\A^Bf — B^A^)j 

i=Wl 


(CS) 


where n is the number of frequency points and the definitions of A^, B^, A^, B^, A^, and B^ are given in equations (8) and (9). 
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